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Abstract 

We study the large A'^ (planar) hmit of pure SU{N) 2+1 dimensional Yang- Mills theory {YM2+1) 
using a gauge-invariant matrix parameterization introduced by Karabali and Nair. This formulation 
crucially relies on the properties of local holomorphic gauge invariant collective fields in the Hamiltonian 
formulation of YM2+1. We show that the spectrum in the planar limit of this theory can be explicitly 
determined in the A'' = 00, low momentum (large 't Hooft coupling) limit, using the technology of the 
Eguchi-Kawai reduction and the existing knowledge concerning the one-matrix model. The dispersion 



relation describing the planar YM2+1 spectrum reads as ui{k) = y + m^, where n = 1,2,... and 
m„ = nmr, where m-r denotes the renormalized mass, the bare mass m being determined by the planar 
't Hooft coupling gyM^ via m = The planar, low momentum limit, also captures the expected 

short and long distance physics of Y M2+\ and gives an interesting new picture of confinement. The 
computation of the spectrum is possible due to a reduction of the Y M2+\ Hamiltonian for the large 't 
Hooft coupling to the singlet sector of an effective one matrix model. The crucial observation is that the 
correct vacuum (the large A*' master field), consistent with the area law and the existence of a mass gap, 
is described by an effective quadratic matrix model, in the large A'^, large 't Hooft coupling limit. 
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1 Introduction 



The study of the large N hmit of Yang-Mills theory is one of the grand problems of theoretical physics. In 
recent years, a new viewpoint has emerged concerning the planar limit of gauge theories, mainly motivated 
by recent insightful advances in string theory based on the duality between string and gauge theories Pp. 
Nevertheless, a precise formulation (a prerequisite for a solution) of the elusive "QCD string" is still lacking. 

In this letter we take a fresh look at this problem in the setting provided by 2 + 1 dimensional Yang Mills 
theory {FM2+1) - a highly non-trivial quantum field theory [2], [3]. This theory is expected on many grounds 
to share the essential features of its 3 -I- 1 dimensional cousin, such as asymptotic freedom and confinement, 
yet is distinguished from its 3 -I- 1 dimensional counterpart by the existence of a dimensionful coupling. We 
regard this study as a stepping stone towards the 3-1-1 dimensional theory. 

Interestingly enough, we are able to make a precise statement concerning the spectrum of this theory in 
the large N, reduced, low momentum limit. In this limit, we argue, the generic features of the Yang-Mills 
vacuum are fully captured by a quadratic large N matrix model. We perform explicit computations in a 
well-defined framework utilizing a momentum expansion of the reduced, planar effective action given in terms 
of the local gauge invariant variables which correspond to the only propagating physical polarization. Our 
approach utilizes many recipes from the large N cookbook (the large N reduction, matrix model technology), 
and yet is seemingly not directly related to the recent advances in the understanding of certain planar gauge 
theories from a string theory (gauge theory/gravity duality) point of view. This of course does not mean 
that a possible indirect connection is non-existent. We note that our approach can be understood as a target 
space, Hamiltonian formulation of an effective string field theory describing the planar 2-1-1 dimensional 
QCD string. 

Our work is crucially based on beautiful results derived in the remarkable work of Karabali and Nair 
[3]. They have provided an explicit gauge invariant reformulation of YM2+1 in terms of local holomorphic 
variables. On a more practical level, Karabali and Nair have been able to compute the string tension in their 
Hamiltonian approach which is in excellent agreement (up to 3%) with the existing lattice data 4 . That 
striking result as well as the computation presented in this letter clearly point out that the Karabali-Nair 
approach has some truly remarkable features which can lead to potentially dramatic results in the arena of 
2-1-1 dimensional gauge theories. 

It is reasonable to believe that if a solution of YM2+1 is to be found at all, it will be in the planar limit. 
Consequently, we consider here the large N limit of YM2+1 using the Karabali-Nair parameterization, and 
this is the distinguishing feature of our work. We consider the spectrum of the theory and are able to show 
that there exists a mass gap set by the 't Hooft coupling. This result extends but is certainly consistent with 
the results of 

The crucial element in the computation of the spectrum is a large N reduction of the YM2+1 Hamiltonian, 
in a well-defined low momentum limit (large 't Hooft coupling), written in terms of the Karabali-Nair 
variables, to the singlet sector of an effective one Hermitian matrix model. The singlet sector is selected 
by the presence of a local holomorphic invariance of the planar vacuum which arises in the Karabali-Nair 
formalism. This reduction procedure enables us to write a self-consistent gap equation for the planar sector 
of YM2+1. What is most important is that our approach provides a well-defined momentum expansion 
of the full effective local gauge invariant collective field theory of YM2+1 in the reduced, planar limit. In 
particular, the correct Yang-Mills vacuum (i.e., the large N master field), consistent with the area law and 
the existence of a mass gap, is captured by an effective quadratic matrix model, in the large iV, large 't 
Hooft coupling limit. This effective theory of gauge invariant holomorphic loop variables can be in principle 
used for other calculations, such as the determination of various correlation functions. 

The outline of this paper is as follows: In Section 2 we review the Karabali-Nair variables and then 
in Section 3 we investigate in detail the Karabali-Nair collective field theory Hamiltonian. The large N 
limit of this Hamiltonian is studied in Section 4 and the planar spectrum of the same in Section 5. A few 
more technical details related to the analysis of the collective field Hamiltonian are collected in a separate 



2 



Appendix at the end of this letter. 



2 The Karabali-Nair variables 

The Karabah-Nair approach can be summarized as follows ^3 : consider an SU{N) Y M2+1 in the Hamiltonian 
gauge ^0 — 0. Write the gauge potentials as Ai — ~it°-A1, for i = 1,2, where are the Hermitian N x N 
matrices in the SU{N) Lie algebra = if°'^H'^ with the normalization 2Tr{t°-t^) — (5"^. Define complex 

coordinates z ~ xi — 1x2 and z — xi + 1x2, and furthermore 2A°- = A° + iA2, 2A'^ — AI — lA"^- 
The Karabali-Nair parameterization is 

A^ -d^MM-^, A = +{M-ydiM'^ (1) 

where M is a general element of SL{N,C). Note that a (time independent) gauge transformation A 
gAg~^ — dg g"^, A — > gAg~^ — dg g~^, where g G SU{N) becomes simply M gM. The variables M 
correspond to holomorphic loops; their most important, and perhaps unexpected, property is locality! The 
corresponding local gauge invariant variables are given in terms of "closed loops" H = M^M. Note that the 
standard Wilson loop operator may be written 

$(C) = TrPexp{-i j dz d^HR-'^} (2) 

and thus is closely related to the local H variables. 

Now one might wonder whether the parameterization Q is well-defined. In fact, the definition of M 
implies a holomorphic invariance 

M{z,z) M{z,z)h\z) (3) 

M\z,z) h{z)M\z,z) (4) 

where h{z) is an arbitrary unimodular complex matrix whose matrix elements are independent of z. Under 
the holomorphic transformation, the gauge invariant variable H transforms homogeneously 

H{z, z) ^ h{z)H{z, z)h\z) (5) 

This is distinct from the original gauge transformation, since it acts as right multiplication rather than left 
and is holomorphic. One way to understand its appearance is that the parameterization ^ can be formally 
inverted in the form 

M{x, x) = ( 1 - / d^z G{x, z)A^{z, z) + ..\ V{x) (6) 



where G is the Green's function, dzG{z,x) = 5^'^\z — x) and V is an arbitrary matrix with only anti- 
holomorphic dependence. The theory written in terms of the gauge invariant H fields will have its own local 
(holomorphic) invariance. The gauge fields, and the Wilson loop variables, know nothing about this extra 
invariance. We will deal with this, as in 3;, by requiring that the wave functions (or cquivalently, physical 
states) be holomorphically invariant. The insistence of the holomorphic invariance of the vacuum is of crucial 
importance for our main argument in what follows. 

One of the most remarkable properties of this parameterization is that the Jacobian relating the measures 
on the space of connections C and on the space of gauge invariant variables H can be explicitly computed 

d^i[C] = o-d^[iJ]e2'=^'^"'^''^[-^l (7) 

where ca is the quadratic Casimir in the adjoint representation of SU{N) (ca — N) and 

Swzw{H) = I d^z TrH-^dHR-^dH + 4r I e^^^TrR-^d^^HH-^duHH^^dxH (8) 

2vr J 12n J 
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is the level —ca hermitian Wess-Zumino-Witten action, which is both gauge and holomorphic invariant, a 
is a constant determinant factor. Thus the inner product may be written as an overlap integral of gauge 
invariant wave functionals with non-trivial measure 

(1|2) = J rf^[i/]e2'=^^"'^»'(^)^'i«'2 (9) 



The standard YM2+1 Hamiltonian 



Tr gi-j,,E.,' + -^B' (10) 
\ Sym J 

can be also explicitly rewritten in terms of gauge invariant variables. The collective field form "S" of this 
Hamiltonian (which we will refer to as the Karabali-Nair Hamiltonian) can be easily appreciated from its 
explicit form in terms of the natural "current" -like gauge invariant variables""^ J = —dzHH~^, 



where 



_ 9ymCA „ , s CA Sab i fabcJ%x) 

ZTT TT^ (a; — y)^ n [x — y) 

The derivation of this Hamiltonian involves carefully regulating certain divergent expressions in a gauge 
invariant manner. We note that the scale m is essentially the 't Hooft coupling. 

At this point we remind the reader about the difference between collective field theory and effective field 
theory. Collective field theory |H] is simply based on a choice of collective variables, appropriate to the physics 
in question. The technical difficulty usually lies in the explicit change of variables, which generically renders 
the collective field theory horribly non-local. (This is, for example, the case with the collective field theory of 
canonical Wilson loop variables.) The crucial requirement that the large N collective field theory has to meet 
is the factorization of vacuum correlators. The factorization in turn implies, by the resolution of the identity, 
that the only state that controls the physics at large TV is the vacuum. Notice that because of factorization 
at large N, one essentially has to be concerned with the appropriate classical phase space of gauge invariant 
observables and their canonically conjugate partners and correspondingly the classical Hamiltonian. The 
expectation values (evaluated using appropriate semiclassical coherent states) of the quantum Hamiltonian, 
lead to the required classical Hamiltonian, which in turn is nothing else but the collective field Hamiltonian 

The truly amazing feature of the Karabali-Nair holomorphic loop variables is that they are local and 
that the corresponding Jacobian can be explicitly computed! This Jacobian, determined in terms of the 
Wess-Zumino-Witten action, enjoys certain analyticity properties which render it unique and independent 
of regularization ambiguities. 

The passage to the collective field Hamiltonian may be thought of as a starting point: having performed 
the change of variables, we may then analyze the theory using effective field theory techniques. In particular, 
one expects in general that there will be renormalizations; because of the local nature of the variables, one 
may expect that a suitable perturbative analysis can be found which sensibly deals with such matters. Indeed, 
we will describe such a formalism here and explain how the dynamics of the mass gap and confinement arises. 

Finally, note that the inner product can be put into a canonical form 



(1|2) = J dfx[H]^l^2 (13) 



provided we perform a redefinition of wavefunctionals $ = e'^-*'^"'^"''^^-'^'. In so doing, there will be a 
corresponding adjustment of the collective Hamiltonian, containing new terms. We will display this explicitly 

'^The J variables transform as connections under the holomorphic transformation. 
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in the following sections, but note here that the most important effect is to add a term of the form^ 
m^Tr(dHdH~^), which will later be understood to correspond to the appearance of the mass gap in the 
large TV limit. Of course, near gvM ^ 0, we would take the perturbative vacuum and this term is of no 
particular relevance — the description of the theory in terms of A and A is adequate. However, this term has 
arisen from the Jacobian of the path integral measure and although the change of variables in the measure is 
essentially an operation on the classical configuration space, we claim that one is lead to address the physical 
non-pcrturbative vacuum. It is with respect to this vacuum that the mass gap appears. 

One notes that this collective field formalism is true for any rank of the gauge group, and in particular 
agrees with the large N 't Hooft counting. We will obtain additional insight into the dynamics of confinement 
however by examining the theory in the large N limit. First, let us continue reviewing the results of Refs. 

m 



2.1 The vacuum wave functional, area law and string tension 

One of the major results of the Karabali-Nair collective field theory approach is the analytic deduction of 
the area law and an explicit computation of the string tension "B . 

The computation of the string tension is achieved by an approximate formula for the vacuum wave 
functional First, one notices that ^' = 1 is annihilated by the kinetic term and is normalizable given the 
non-trivial measure (due to the normalization of the WZW path integral) . 

One may find a vacuum wave functional annihilated by the total collective Hamiltonian (Ti^' = (T + 
V)'^ = Q) by expanding this equation in powers of (roughly) B /m? . To leading order, one finds 



~ exp 



^ ^ B{x)(^^==]B{y) 



(14) 



Note that this wavefunctional apparently interpolates between the low and high momenta regions. At high 



momenta, this wavefunction correctly has a form corresponding to free gluons 'I' e ^"ym ' ' ^ appro- 
priate to the conformally invariant two-point function of gluons, (AA) ^ gYAi/lM- the low momentum 
region, the momentum factor is cut-off, and B'^/k ^ B"^ /m. Although higher order corrections are expected 
to be non-local, this can make sense self-consistently, if the theory may be re-organized into an expansion in 
inverse powers of m. 

As explained in the low momentum limit 

provides a probability measure vl/*v[/ equivalent to the partition function of the Euclidean two-dimensional 
Yang-Mills theory with an effective Yang-Mills coupling g^j-, = mgYM- Using the results from 0, Karabali, 
Kim and Nair deduced the area law for the expectation value of the Wilson loop operator 

($) - exp(-o-A) (16) 

with the string tension following from the results of [H] 

4 iV^ - 1 

= 9ym (17) 

This formula agrees beautifully with extensive lattice simulations 4 , and is certainly consistent with the 
appearance of a mass gap. Notice that this result is once again in full agreement with the large N 't Hooft 
expansion. 



This comes from a piece of Swzw ■ 
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3 The collective field Hamiltonian 



We are interested in studying the planar limit of the Karabali-Nair approach to YM2+i- The large N limit 
is expected to be controlled by a constant cx3 x cx3 matrix configuration called the "master field" [7] . Such a 
configuration should capture correctly both the short and long distance properties, that is, both asymptotic 
freedom as well as confinement. As already mentioned above, knowing the master field configuration is, in 
a very precise sense, equivalent to knowing the correct vacuum at large N, which is the most remarkable 
result of the Karabali-Nair approach as seen from the preceding section. Small perturbations around this 
configuration should lead to the spectrum of glueballs. In the planar limit this spectrum is expected to be 
equidistant, consisting of an infinite number of non-interacting massive colorless particles. 

Can we compute the planar spectrum of YM2+1 using the Karabali-Nair scheme? The claim of this letter 
is that in the low momentum, or equivalently large 't Hooft coupling limit, the planar spectrum of YM2+1 
is explicitly computable. Because the knowledge of the master field is equivalent to the knowledge of the 
true vacuum in the large N limit, the crucial property to be used is of the holomorphic invariance of the 
gauge invariant variables H — s- VHV . This fact, in combination with the known properties of the spectrum 
of singlet states in the one matrix model is what makes the computation of the planar spectrum possible. 

The crucial observation we make here is that the above Karabli-Nair vacuum wave functional consistent 
with the area law (and as we will see, with the existence of a mass gap), is captured by an effective quadratic 
matrix model, in the large N , large 't Hooft coupling limit. The usual power counting arguments (based on 
the power expansion in terms of the large N 't Hooft coupling) can be applied to the Karabali-Nair vacuum 
wave functional by assuming a WKB ansatz, — exp(r), which is consistent with factorization at large N. 
By assuming a local expansion of V in terms of gauge invariant observables, such as J currents, (involving 
possible non-local, J-independent kernels), the leading term in the large N 't Hooft coupling is given by a 
quadratic expression in terms of currents! We will show that this Karabali-Nair vacuum wave functional is 
reproduced by an effective quadratic matrix model involving renormalized couplings (such as the mass m). 

In order to get to this result our first aim is to better understand the structure of the Karabali-Nair 
collective field Hamiltonian. That is the subject of the present section. In the following section we will study 
the collective field Hamiltonian in the planar limit. 

The Karabali-Nair collective field Hamiltonian can be thought of as a string field theory Hamiltonian for 
a pure QCD string in 2 + 1 dimensions. Indeed, the Karabali-Nair variables represent local gauge invariant 
variables and act operatorially on a true, non-perturbative Fock space of 2-1-1 dimensional Yang-Mills theory. 
Thus this description is intrinsically second-quantized and gauge invariant, and as such does qualify as a 
QCD string field theory. Note that this second quantized theory does act as an interacting theory in the usual 
space, in other words this QCD string field theory is not formulated on a loop space, precisely because the 
string field in this approach (identified with the Karabali-Nair variables) is local. As we will see in the next 
section, the Hamiltonian of this effective QCD string field theory is generically non-local. A first quantized 
worldsheet theory remains for the moment elusive; presumably such a theory is interacting, and moreover, 
can be deduced from the target space second quantized theory studied in this letter. Nevertheless, a purely 
first quantized description should not be considered satisfactory as it would yield just the spectrum and 
compute S-matrix-type observables. The second quantized formulation in principle contains full knowledge 
of the non-perturbative Fock space. In this letter we concentrate on the form of the non-perturbative planar 
vacuum and the spectrum of gauge invariant excitations around it. 

3.1 Hamiltonian 

The classical mechanics of the H variables is somewhat complicated by the constraint det H — 1. The 
Hamiltonian is more easily found using traceless variables - for example, it is convenient to use the currents 
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J = £AdHH-^. The Karabali-Nair Hamiltoniangj is 
where 

n^\x,y)= (19) 

7r^(x — yj^ Tr[x — y) 

This is derived within a consistent gauge-invariant regularization scheme. The last term in (|18|) is the 
potential term, and follows from the precise relation dJ — BAI~^ . 

Alternatively, suppose we consider expanding around the constant solution as 

H = e'^, (20) 

where ip — ip"'t°' is Hermitian and traceless. The perturbative vacuum is described by iJ = 1, and we are to 
expand the theory in powers of Lp. Given the parameterization of H , one finds 

H-^dH = d^ + \ [d^, ^] + i [[d^, (p] + . (21) 

Using the (adjoint) notation'^ ip"-^ = ip'^f"-^'^ (or equivalently, y'"=(^'>c ~ CAf'^), we find 

H-^dH ed(p''eba[(p] (22) 

where'* ef,a[</'] = ^ba — ^^ba — ^{f'^)ba . . • is a functional of ip. Note that in this notation, we have a generalized 
non-Abelian bosonization formula 

= !^eafc[^]a^^ (23) 

Moreover one can also establish the following useful formula (see the Appendix): 



SJ-{x) CAjy "^(^■^)<5(^^(y) 

where D'ab ^ Sabd - jfabc^^d + ifabcdip" + . . .. 

Given this dictionary between J and (p variables, the collective field theory Hamiltonian may also be 
expressed in terms of the ip> variables. It has the form 

-HknI^] = [p'^Mx)-4^+f Q^'Mx,y) I I +^ [ gj°Mgj°M (25) 
Jx bip°-{x) J^ y dip'^ix) dip>'{y) mcA Jx 

The formulae for the fmictionals P°[</'] and can be found in the Appendix. Formally, we find point- 

split versions of: 

P^Wx) = j{e-\x)UG(,^,)H{zUG(,^^){e-\x))da,c (26) 

g'^^M(x,y) = _^^(e-i(y)),,G(,,,)i?(z)edG(,,,)(e-i(x)),, (27) 

As in the Appendix, H^b is the adjoint representation (e~"^)af,. As we will see, the change of variables from J 
to </? does introduce extra technical problems in the regularizaton of various expressions, but these issues do 
not effect the final physical result. In particular, the mass gap is not an artefact of regularization, precisely 
because the collective Hamiltonian was derived using a consistent gauge-invariant regularization scheme. 

^We distinguish the two representations by exphcit indices in the following formulae when necessary. Note that here 1/3°'' is 
real skew-symmetric. 

''There is an explicit resummation ej^j = it/j"-*^ (e"'*' — 1), in the adjoint notation. 



7 



3.2 WZW 

Next, let us look at the WZW action. We have 

Swzw = -^l Tr H-^dH.H-^dH + I ef^^^Tr H-^d^H.H-^d^H.H-^dxH (28) 
2-K Jm 127r Jp 

where dP = M and H is Hermitian. Given the parameterization of H, we can then evaluate 

Tr H-^dH.H~^dH = ^dip''d^''eac[v]ebc['p] (29) 
= ^d^'^B^'g^M (30) 

where g[ip] = e[<^]e[_^] = / - ^if"^ + ... 
The WZW term may be written 

J dhds /la6cMa(p"Vc^.¥'' (31) 

where hate = fdefSadete^cf- This term is of course a total derivative, and gives 

^ J dh dif^d^Xaiif] (32) 

where habc = bab,c - bcb,a - bacb- {b = \ip + . . . IS antisymmetric). 

Thus we arrive at the familiar result that the full WZW action can be written as a sigma model action 

Swzw = / d^'^d^'GZ^'^iv) (33) 

where 

GZ^'^{v) = {g + ^b)M (34) 

For small fields, which correspond to H of order one, one can use the standard background field method. 
We obtain 

GZ'^{^)={^ + l^-Y2^' + --) (35) 



ba 



The calculations presented below can be obtained by following this formalism. 

Now, we may proceed to do the redefinition of the wavefunctionals. The effect of this redefinition is to 
give a new collective Hamiltonian, which is a similarity transform of Hkn 



where 



and 



n' = e^'^^wzw'We-^'^'^"'^^ = W2 + Wi + Wo (36) 

W2=/ Q'''Mx,y) ^ ,^ . (37) 
Jx,v Sifi'^ix) dip^y) 



P'^M(X)-2C^ / Q«^[^](^,y)^g: 



(38) 



and, finally 



Wo = -— f d{eabd^'')dieacd^n - CA f P"M(^, . , 
-CA I Q^'W^^y) 

Jx,y 



x)^-p^ (39) 



CA 



Sifi°'{x)Stp''{y) S(p°'{x) 5ifi^{y) 



(40) 
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These are formal manipulations and one needs to take care to regulate the expressions in a gauge and 
holomorphic invariant manner. In fact, in general this expression will be non-local; however, as we will see 
later, the real utility of the ip variables is that the Jacobian discussed above essentially generates a mass gap. 
The mass gap is set by the 't Hooft coupling itself, and so there should be a self-consistent low-momentum 
expansion, in powers of k/m. 

Let us work in an expansion in powers of (p. Recall that eba — ^ba — ^fba — ^(<P^)&a • ■ • and Gab = 



(g + ib)ab = Sab 



SSwzw 



Note that we have 



1 
1 

2^ 



ddip'^ + -fabcdif^Bif^ 



+ Gca.dM - Gcd,aVp\) 



Using this expression, we may derive, as an expansion in (p, 



9ym 



25 



YM 



9ym 



dip''i-dd)dip'' 



(41) 



(42) 



where C{x,y) is in general a non-local kernel whose inverse is, to leading order in k, 

This collective Hamiltonian TC' is also in general non-local. The terms given in eq. (|42|) represent the part 
of the Hamiltonian which controls the ground-state in the planar limit; we note however that the parameter 
m should be understood as a renormalizcd coupling to^- The quadratic terms determine a dispersion relation 
of the form 



A(fc) = |fc|2(£;2_|fc|2_ 



(43) 



This indicates that the field (p is not canonically normalized. However, since (p has no physical zero-mode^, 
the transformation to canonically normalized excitations is non-singular, albeit non-local. We will explore 
this further in the next section. 



4 N = oo reduction 

All of the above analysis is valid at any N. Let us now consider the theory in the planar limit. For self- 
consistency of our presentation we briefly review the Eguchi-Kawai reduction [S], which is strictly valid in 
the N = oo limit. We discuss the reduction in the continuum and concentrate on the matrix scalar field 
theory, obviously relevant for our discussion. 

Consider the following general local action for a scalar matrix field (p 

Siip) =nJtt [^ida^r + Viipix))] (44) 

The N = oo reduction is captured by the following recipe: 

1) Replace the position dependent scalar matrix field (p{x) by 

(^(x) ^ e^^-^V^e-'^"^" (45) 

and 

tt{x) -> e'-P-^Vfle-*^"^" (46) 
where Pa is a diagonal Hermitian matrix P"^ = diag{pi,p2...p%) and 

2) replace the derivative operation da^{x) by 

daVix) -^i[Pa,V^R] (47) 



^Constant ip is equivalent, via a constant holomorphic transformation, to i/p = 0. 
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and 

3) finally, replace the continuum action per unit space-time volume by 

SR^Tr[~]^[Pa,^RY + VR{^R)] (48) 

Then the correlation functions of the continuum theory (in the N ~ oo limit) are given by the correlation 
functions of the reduced theory, after integrations over the eigenvalues of the momenta pi 

N 

{FMx)]) ^ I l[d^MF[e^'''^^^Re'^P^^^]R) (49) 

i—l 

Following this prescription we can easily write a reduced form of the local part of the collective field Hamil- 
tonian in terms of (p variables. This is sufficient, because this part of the collective field Hamiltonian controls 
the planar vacuum at large 'tHooft coupling, as we will see in what follows. 

As we would like to go to large N , we consider the continuum Eguchi-Kawai reduction of the approximate 
local expression of the Hamiltonian H42() . The low momentum limit 

% « 1 (50) 

defines the large 't Hooft limit, given the fact that the gauge coupling in 2 + 1 dimensions is dimensionful. 
As we noted above, in the large iV, large 't Hooft coupling limit, the vacuum wave functional is given by 
an exponent of a quadratic functional in the current variables. We will show that the large N large 't 
Hooft coupling limit is self-consistent in the sense that it leads to the correct vacuum implying a gap in the 
spectrum. 

To proceed, we introduce the following change of variables in momentum space (a non-local change of 
variables in real space) that we alluded to in the previous section, in order to get a canonically normalized^ 
kinetic term 

r{k) = J^^^{k). (51) 

V 9ym 

We see that in the low momentum limit of the reduced Hamiltonian we simply get 

Tr (^2 + „^2^2 _ ^[p^^ [pa^ ^ _ (52) 

where the canonical momentum tt — 

As it stands this Hamiltonian apparently describes N"^ — 1 massive degrees of freedom with a mass 
proportional to the square of the gauge coupling, as noted originally in |2] . Obviously this does not seem to 
give the confining spectrum we expect of YMi^\\ As a matter of fact, in the limit of the zero Yang-Mills 
coupling the familiar perturbative spectrum of gluons is readily recovered, as we discussed in Section 2. 
This result is natural as we have expanded Hix) — exp((/?(a;)) around the "perturbative vacuum" H — \ to 
quadratic order in The new interesting feature in this discussion is the presence of the gauge invariant 
mass term, whose origin was the Jacobian of the transformation to these variables. 

Nevertheless we claim that one can gain important insight into the nature of the planar vacuum provided 
we remember that the master field, in terms of gauge invariant H — exp((^) variables, is supposed to 
be a constant c» x c» matrix that transforms homogeneously under the residual constant transformations 
H — !■ hHh^. In particular, the vacuum is preserved by constant unitary transformations. Therefore in the 
planar limit the matrix cj) will also be an cx3 x cx3 matrix that transforms as cj) h(f)W . We must require 
invariance under constant unitary transformations on the zero-mode of this field.'' 

®Notc that wc normalize to 1 rather than A'^ here for simphcity. 

^This is a somewhat subtle point beeause of the momentum-dependent change of variables I51i . 
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Thus, the dynamics of the planar master field of YM2+1 which describes the vacuum in the planar, low 
momentum limit, in the well-defined momentum expansion of the reduced Hamiltonian, given in terms of 
matrix variables, is determined by the following one matrix model Hamiltonian 

i 1 Tr (tt^ + ml4>^ + ...) ^ Tr (^-^ + ml4>^ + (53) 

The vacuum in this approximation is captured by the singlets, invariant under the residual unitary transfor- 
mation (j> —> h(f)h'^ . This is a huge reduction in the number of states in the large N limit, described by the 
density of -/V — *■ 00 eigenvalues of the matrix ^! 

Therefore, due to the holomorphic invariance of the vacuum, the spectrum of the planar limit of Y M2+1 
is determined in the planar, low momentum (large 't Hooft coupling) limit, by the spectrum of the singlet 
states of this effective one matrix model. 

Before proceeding further, let us note that in order to compute correlation functions we would need 
to include momentum dependent terms in the reduced Hamiltonian in terms of matrix variables. One 
obvious technical complication one has to deal with in order to get the leading momentum expressions for 
various correlation functions is that the momentum-dependent terms in the reduced Hamiltonian are not of 
the usual kind considered in the literature of the one matrix model (i.e. they are not traces or multi-trace 
terms involving a single matrix) . Consequently, the computation of correlation functions will be significantly 
harder, although we believe the formalism presented here is sufficient for a discussion of the spectrum. In 
this letter we restrict our attention to the spectrum. 



4.1 The vacuum wave functional, one more time 

Now we demonstrate that the effective quadratic matrix model captures the correct physics of the vacuum. 
First, we can easily find the correct vacuum wave functional (which leads to the area law and a successful 
empirical expression for the string tension) using the above matrix field theory approach. Given the quadratic 
matrix field Hamiltonian H52|l. we see that the ground state wave functional is a Gaussian 




To compare to the previous discussion, we should consider the effect of transforming back to the non-trivial 
inner product of the original collective field theory. This amounts to a simple shift of the exponent of the 
wave functional 

* = exp ^-i j (j){x) [~-mr + ^ml - </>(j/)^ (55) 

Conversion to the J variables to linear order in </> fields gives a wave functional which perfectly matches the 
expression for the wave functional l|14|l . which, as we have seen above, leads to successful predictions of the 
area law and the string tension. 

Thus we explicitly see that in the large A^, low momentum, that is, large 't Hooft coupling limit, the 
effective quadratic matrix model leads to the correct physics of the Yang-Mills vacuum, as described by a 
vacuum wave functional consistent with the area law. This result should not come as too much of a surprise, 
since it is essentially the same calculation that we considered in Section 2.1. However, it is a useful check 
that the reduction has not thrown away anything important about the vacuum. 

We are now ready to show that the same effective matrix model describes the correct spectrum of 
excitations about the vacuum of the large 2 + 1-dimensional Yang-Mills theory. 



5 The planar spectrum of YM^^i 

In this section we determine the planar spectrum of YM2+1 from the spectrum of the singlet states of a 
one-matrix model. The one matrix model is a well-studied system [5]. It can be understood from many 
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points of view of the known large N technology |S] 1^1201 • For example, it is well known that the planar limit 
of a quadratic matrix model in the singlet sector is completely captured by the semicircular Wigner-Dyson 
distribution of the density of eigenvalues (A) of the matrix (j) 

PwdW - -y/2a-mf.\^ (56) 

TT 

where A G [—A, A], A being the point where the square root vanishes, that is A = (a is just a Lagrange 
multiplier associated with the normalization condition satisfied by the eigenvalue density. In the equivalent 
fermionic formulation, it is the Fermi energy.) Furthermore, the master field of a Gaussian matrix model I 10| 
can be described in terms of noncommutative probability theory II 1| (for which the semicircle distributions 
play the role completely analogous to the Gaussian distributions of the commutative probability theory) as 
an operator (j) = a + acting on a Fock space build out of a and operators and the vacuum |0), a|0) — 0. 
a and satisfy the Cuntz algebra aa^ — 1. 

The spectrum on singlet excitations can be easily determined by perturbing around the ground state 
value, determined by the semicircle law, assuming e"^* time dependence and reading off the normal modes 
from the resulting wave equation. (For example, this has been done in detail in for the case of the 
one-matrix model.) Because these calculations are crucial for our main claim, we review them from a couple 
of different points of view, mainly following the work of jl2| . 



5.1 Spectrum of singlets from the collective field theory 

As stated above, the density of eigenvalues p(A) is the correct variable that describes the dynamics of the 
large N one matrix model. One introduces first pk = Tr eKp(—ik(j)) , and then defines the collective field p(A) 
as the Fourier transform of pk- Here p(A) is positive definite and normalized as J p{X)dX = N. 

Suppose that the one matrix model is described by the general Hamiltonian Tr{^Tr^ + V{(j))). The 
collective field Hamiltonian for p(A) and its conjugate momentum H(A) reads as follows 

i J dA[p(A)(9AH(A))2 + + Vi\)p{X) - ap[X)\ + aN (57) 

The classical equations of motion for the collective field theory are 

p(A) - -aA[p(A)aAn(A)], ri(A) = \{dm\)f + \i^p{\f + V{\) - a (58) 
The solution of these equations that satisfies the positive definiteness of the eigenvalue density is such that 

9An(A)-0, p{\)^pwd{.X)- (59) 

The spectrum can be now determined by small perturbations around this classical ground-state solution. 
In particular, introduce following ^21; = d\Tl{X) and p(A)^ = Pq{XY + and expand the Hamil- 

ton equations (|58|l to first order in ^ and rj. Assuming the periodic time dependence exp(— iwt) for both 
perturbations, we may eliminate ^ to get the following wave equation for 77 



+ || = (60) 



where q{X) = j^[2{a - V{x))]-^2dx. 

The boundary conditions are determined from the constraint on the density of eigenvalues / p{X)dX = N 
and the vanishing of po for A = ±A, from which one infers that §| = for q — q{A). This implies that 

the wave equation 1)60(1 describes a finite length string with free ends (for example, rj ^ cos ("-""^(aj)) fro™ 
which we infer the spectrum 

LU = nujc (61) 
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where the fundamental frequency is given by 




-V{X)r-^j . (62) 
In the case of a quadratic potential V{x) — , this can be integrated precisely into a gap-like equation 



nir (63) 



Of course, the same result follows by directly perturbing the collective field Hamiltonian and determining 
the normal modes of the quadratic part [12] . 

As is well know, the above result (|61|l can be readily understood from a different, and somewhat more 
intuitive, point of view. As is well known [H], 5S| the singlet sector of the one matrix model can be understood 
in terms of an effective system describing N free spinless fermions in a quadratic potential. The relevant 
Hamiltonian is 

+ V^(A) (64) 

The ground state in this picture corresponds to a filled Fermi sea, with the corresponding Fermi energy a 
so that 

N = - [ dX{2{a^V{X)) (65) 



TT _ 

The density of states in the semiclassical limit is given by 

^ J dXdpSia - (ip2 + v{X))) = dX{2{a - t/(A))-^ (66) 

For a quadratic potential ^(A) = ^ A^ the fundamental frequency ujc — mis given as an inverse of the density 
of states at the Fermi level. Of course, this result agrees with the corresponding collective field theory result. 
Notice, that lOc is nothing but the classical frequency of the corresponding classical trajectory describing the 
motion of free fermions in the effective potential V{X). Finally, the level degeneracy is also easy to discuss 
in the fermionic picture. 

The crucial point here is that the singlet spectrum of a one matrix model has a gap, determined by the 
mass parameter m, is equidistant and consists of free non-interacting excitations. Going back to the planar 
limit oiY M2+1, we see that for the large value of the 't Hooft coupling, the master field of the planar YM2+i, 
a constant 00 x 00 matrix left invariant under the residual holomorphic transformations, is equivalent to the 
eigenvalue distribution of the singlet sector of a quadratic one matrix model. Therefore the spectrum has a 
gap determined by the mass m and is equidistant. The mass gap emerges because of the finiteness of the cut 
of the semi-circle distribution and the nature of the boundary conditions at the end points of the cut. This 
offers an interesting view on confinement. 

The Karabali-Nair formalism is completely consistent with Lorentz invariance as shown in ^j. The 
Lorentz invariant form of the planar YM2+1 spectrum is then given by the following dispersion relation 
describing an equidistant spectrum with an explicit mass gap 



uj{k) = ^^ + ml (67) 

where n = 1, 2, to„ = nm^, where is the renormalized gauge invariant mass and the bare mass m^ is 
determined by the planar 't Hooft coupling gy^^ ™ = '^^iif ■ First of all, this formula elucidates the 
meaning of the parameter m in the planar limit of YM2+1 in the large 't Hooft coupling limit. Note that 
as the 't Hooft coupling is taken to zero we recover the dispersion relations for massless particles. This is 
in accordance with the recovery of a Coulomb potential at short distances, as shown in 3 . Finally, as we 
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have seen, the matrix model knows about the vacuum wave functional which in the high momentum limit 
does describe a theory with a Coulombic potential and in the low momentum limit leads to the area law for 
the expectation value of the Wilson loop and an explicit formula for the string tension which matches the 
numerical data. 

5.2 Matrix model reduction and the planar spectrum 

Obviously, the structure of higher order momentum dependent terms in the effective gauge invariant Hamil- 
tonian is in general non-local and very involved. The question that we would like to briefly consider is how 
the higher order terms influence the leading quadratic result for the mass gap and the spectrum. Even 
though the higher order corrections are explicitly non-local, we will argue that the existence of a mass gap 
in the large N, large 'tHooft coupling limit, implies a self-consistent expression for the mass term even upon 
the inclusion of the higher order terms in tp. 

As we have seen, the planar vacuum at leading order in the 'tHooft coupling is governed by the quadratic 
part of the collective field Hamiltonian, and the higher order contributions to the wave functionals can be 
re-organized into an expansion in inverse powers of m. Motivated by this observation, suppose we cut-off the 
momentum integrals in all non-local terms in the collective field Hamiltonian by the mass gap M and then 
expand in momentum and finally perform the large N reduction. This procedure is difficult to implement 
technically, yet nevertheless, the general form of the effective reduced Hamiltonian involving (j) matrices, 
would look as follows 

+ + (68) 

The coefficients in V{<j)) would only involve powers of the renormalized Karabali-Nair mass parameter m, 
denoted by rrir- 

Now, by insisting on the holomorphic invariance of the vacuum, as in the quadratic case, we see that only 
the singlet sector of this general one matrix model, corresponds to the vacuum of YM2+i. The spectrum 
of the singlets has been discussed in the previous section. The fundamental frequency, corresponding to the 
gap of YM2+1 is given by 

A(m^) 

[2{a~V{\mr))]-^d\)-^ (69) 

By choosing Wc to correspond to the physical mass gap M, and by evaluating the renormalized mass at the 
scale determined by M, we get a self-consistent gap equation 

^A(m,(M)) 

M=-(y^ [2(a-y(A,™,(M)))]-^dA)-i (70) 

Because of the positive definite nature of the potential, a solution to this equation should exist. The lowest 
numerical value, in the units of the 't Hooft coupling should correspond to the physical gap of FM2+1. 

One important point here is that the lowest order, quadratic contribution is not misleading when we try 
to capture the correct long distance and short distance physics in the planar limit. The ultimate reason for 
this is that the gauge invariant H variables are local, and that the inner product on the space of states is 
computable in terms of these variables. The lowest, quadratic order for the master field (described by the 
dynamics of an effective quadratic matrix model) obviously captures an important piece of the correct ground 
state (i.e. the master field). The effective matrix model leads to the ground state that is also consistent with 
the form of the vacuum wave functional obtained in [3], predicting the string tension which turns out to be 
in excellent agreement with the available lattice data In our case, the numerical value of the mass gap 
as given by the spectrum of singlets of the effective matrix model has to be fitted to the lattice data. 
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5.3 Towards the QCD String 

Finally we add a couple of comments regarding the relevance of our discussion for a possible string theory 
of the 2+1 dimensional Yang Mills theory. 

Perhaps the two obvious questions regarding our results in view of the usual intuition connected to the 
QCD string are ®: 1) how would one demonstrate, within our framework, the expected Hagedorn behavior 
of the density of states and also 2) how would one establish the existence of Regge trajectories expected 
from a QCD string? 

To approach answering these questions recall that the Karabali-Nair variables capture all the degrees 
of freedom of the 2 + 1 dimensional Yang Mills theory. Also, the Karabali-Nair Hamiltonian is the exact 
collective field Hamiltonian. Thus, in the confined phase of the theory one expects the usual 0(1) scaling 
of the free energy and in the deconfined phase the 0{N'^) behavior. Indeed, the effective matrix model we 
have used to demonstrate the existence of the mass gap, corresponds to the counting of degrees of freedom 
expected from a confining phase. On the other hand, the original Karabali-Nair collective field Hamiltonian 
does describe 0{N^) perturbative degrees of freedom. The presence of an explicit mass gap, indicates 
that the self-consistent effective scalar matrix field theory derived from the Karabali-Nair collective field 
description is cut off at the scale determined by the gap M, which would be in accordance with the expected 
confinement /deconfinement transition. 

In order to really establish the stringy Hagedorn behavior (and also demonstrate the existence of Regge 
trajectories) in the planar YM2+1 one needs to attach the Lorentz indices to the oscillators of the singlet 
sector of the effective matrix model. The effective matrix model already describes the correct vacuum in the 
planar limit and the oscillators acting in the corresponding Fock space should obey the algebra of the usual 
oscillators of the free string theory, because of the large N factorization. This we think is the key to deriving 
a QCD string field theory for the 2+1 dimensional Yang Mills theory and the associated stringy (such as 
Hagedorn and Regge) features. Work on this important issue is in progress. 

6 Concluding remarks 

To conclude, the spectrum of physical excitations in the planar limit of Y M2+1 can be deduced using the 
remarkable local holomorphic variables of Karabali and Nair taken in conjuction with some well known 
results from large N master field technology. The analytic understanding of the spectrum is possible due 
to a reduction of the YM2+1 Hamiltonian for the large 't Hooft coupling (low momentum limit) to the 
singlet sector of an effective one matrix model. The huge reduction of the degrees of freedom to the singlet 
sector of the one matrix model is a consequence of the holomorphic invariancc of the YM2+1 vacuum in 
the Karabali-Nair representation. Note also that the matrix model captures the form of the vacuum wave 
functional which leads to the area law and a successful empirical expression for the string tension, reinforcing 
the self-consistency of our approach. 

Obviously the approach presented in this letter is just a first step in the direction of unraveling the 
full planar limit of 2 + 1 dimensional Yang-Mills theory. Given our result we believe that we are now well 
motivated to study the full matrix collective field theory for the gauge invariant holomorphic loop variables. 
We also believe that some of the old as well as currently pursued ideas pertaining to the subject of the 
QCD string jj^l such as the question of integrability, the analogy between gauge and chiral fields, as well 
as target space understanding of the gauge theory/gravity duality should be re-examined in the context of 
2 + 1 dimensional Yang-Mills theory. The gauge invariant collective field Hamiltonian can be understood 
as a string field theory of YM2+i. It would be obviously very interesting to understand the world-sheet 
structure underlying this target space description in order to get even closer to the (perhaps not so) elusive 
QCD string. Last, but not least, we believe that a new way is open for a rational analytic approach to the 

*We thank Ofer Aharony, David Berenstein and Shiraz Minwalla for discussions regarding these issues. 
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large N 3 + 1 dimensional Yang- Mills theory. 
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Appendix 

Here we collect some notation and results. First, define G(^x,v) and G{x,v) 

= 5^^\x - y), BG^x^y) = 6^^\x - y) (71) 

= ^(^'^'^^^ ^^^^ 

Gauge covariant versions are 

^("4) = (9 + ^)r4) = M(a;)M-i(y)G(,,,) (73) 

^("4) = id + ^)r4) = M-t(x)Mt(y)G(,,,) (74) 

We thus find 

SM{x) = - [ D-]^^5A{y)M{y) = -M{x) [ G^,^y^{M-'5AM){y) (75) 

5M\x) = [ M\y)5A{y)D-l^^ = f {MUAM-^){y)G(y^,^M\x) (76) 
Jy Jy 

SJ = -£^(MtMM-t)(a;) + ^ / (MtMM-t)(2/)a,G(,,,) - / [J(x), (MtMM-t)(y)] G(,,,) (77) 
^ ^ ^ y ^ y 

This equation is used to derive Hkn[J] together with the regulated expression Tr t°'D~^{x,x) = ^Tr t"'{A — 
M-'idM''). 

It also follows that 

SH{x) = -H{x) j G(,^y){M-HAM){y) + f {MHAM-^){y)G(y^,)H{x) (78) 
Jy Jy 

But we also have 

5H = HeSip\e[^])ba = eH5^\e^^})ab (79) 
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and so 

Sip'' = 2(ef-^|)baTr t'H-'SH - 2{e^^^)abTr t''SHH-^ (80) 

Thus we conclude 

= 2i(ef;|),,(x)G(,,,)Tr rM-i(y)i''M(2/) ^ -zA/''=(y)G(,,,) (ej-^),,!^) (81) 
= -2^(ef-4),,(x)Tr t^Mt(y)t''M-t(y)G(,,.) ^ *(ef4),e(x)G(.,,)Mt^''(y) (82) 

M^(J)S^\z) ^ *^(y)''G(,,.)(e-i).,,(e^i),eG(.,)Mt(z)^^ (83) 

Introducing a gauge-invariant point-splitting procedure (insertion of a point-splitting Wilson line), we can 
use these expressions in the evaluation of the Hamiltonian 



5yM 



M=(z)M^(x) i5(/9"(x) M'^(z) (5^'=(x) (5<^''(w)(5(/j''(j/) 



<5 /■ „, . 5 



2 



(84) 



(85) 



As indicated in the text, another way to proceed is to work with TLknIJ] ^'^'^ convert to the variables. 
We then need 

= — [eacd ~ ieab,c - eac,b)dip''] -4r (86) 
dip'' ITT ■' dJ°- 

ITT Oj°- 

(note that D' ah ^ ^abd - ^fabc^P^d + ifabcdip" + ■■■), so 



SJ-{x) CAjy ''^(^■^)<5^-(y) 
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